Relaying phase synchrony in chaotic oscillator chains 
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We study the manner in which the effect of an external drive is transmitted through mutually 
coupled response systems by examining the phase synchrony between the drive and the response. 
Two different coupling schemes are used. Homogeneous couplings are via the same variables, while 
heterogeneous couplings are through different variables. With the latter scenario, synchronization 
regimes are truncated with increasing number of mutually coupled oscillators, in contrast to ho- 
mogeneous coupling schemes. Our results are illustrated for systems of coupled chaotic Rossler 
oscillators. 
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I. INTRODUCTION 

The many flavours of synchrony that emerge as a con- 
sequence of different coupling scenarios have been exam- 
ined in detail over the last couple of decades. The roles 
of topology and nonlinearity have been studied, and a 
fair understanding of the different forms of synchrony 
that can arise — given a specific scheme through which 
dynamical systems interact with each other — is under- 
stood to some degree [![. Synchronization in various 
forms is common in forced and coupled nonlinear sys- 
tems 043- The manner in which the interacting sub- 
systems are coupled plays a crucial role in determining 
which form of synchrony arises. For example, when iden- 
tical nonlinear systems are coupled unidirectionally with 
one subsystem (the master) driving the response subsys- 
tem (the slave) complete synchronization occurs. When 
the coupled systems are not identical, generalized syn- 
chronization can occur @, Q . "Mixed" synchrony is ob- 
served in counter-rotating coupled oscillators [9| while 
phase synchronization occurs in mutually coupled chaotic 
systems: the phases of interacting systems are entrained 
while the amplitudes remain uncorrelated @. This is a 
phenomenon of great interest due to potential applica- 
tions in different fields, ranging from physics, chemistry 
to biological and medical sciences ficj . 

Since the coupling can be uni- or bi-directional and 
can be linear or nonlinear II If . and can be through sim- 
ilar or dissimilar variables |l2l . Il3j , as the number of in- 
teracting components increases, the possible variations 
grow exponentially. Our interest in present work is to 
examine a small number of "typical" patterns or motifs 
of coupling, and investigate the different patterns of syn- 
chronization phenomena that result. An additional mo- 
tivation arises from the fact that in a variety of natural 
systems that are subject to forcing, the modulation can 
be either direct, namely when a given system is itself sub- 
ject to driving, or indirect, when it is coupled to another 



system which is the one that is being modulated. Such 
indirect modulation is likely to be operative in biologi- 
cal phenomena fl4j or in networks of coupled dynamical 
systems [l5j . 

In the present paper we study coupled chaotic os- 
cillators which are externally forced. We find that as 
the number of mutually coupled oscillators increase, the 
phase synchronization (PS) regime gets truncated if the 
coupling is heterogeneous, namely through different vari- 
ables, while the same does not hold if the coupling is 
homogeneous, namely through the same variables. 

The different coupling schemes are discussed in the fol- 
lowing Section II where we also study phase synchroniza- 
tion between the drive and the response. The effect of the 
drive is further examined in Section III. The measure we 
used to determine phase synchronization is based upon 
the variation in phase difference with time: this is dis- 
cussed in an appendix that follows the concluding Section 
IV which presents a discussion and summary. 
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II. THE COUPLING PATTERNS 



We first consider the model system of three oscillators 
coupled as schematically shown in Fig. 1. Oscillators 
O2 and O3 (denoted by the subscripts in the variables) 
are diffusively and mutually coupled, and are driven by 
oscillator 0\. The response oscillators O2 and O3 are 
identical and are distinct from oscillator 0\, there is a 
parameter mismatch in the frequencies. The scheme in 
Fig. 1(a) is termed "heterogeneous" since the driving is 
effected through a variable that is not involved in the cou- 
pling, while that in Fig. 1(b) is termed "homogeneous" 
since the driven and coupled variables are the same. 
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FIG. 1: (a) Heterogeneous coupling: the (nonidentical) 
chaotic oscillator Oi is coupled to two mutually coupled iden- 
tical oscillators through one of the variables, here y, while the 
identical oscillators O2 and O3 are symmetrically and bidi- 
rectionally coupled through the x variable, (b) Homogeneous 
coupling: all the interactions are through the same variable, 
here y. 
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FIG. 2: Schematic phase diagram for the phase synchroniza- 
tion states as a function of ei and £2 between (a) oscillator Oi 
and the directly forced response oscillator O2, (b) 0\ and the 
indirectly driven oscillator O3, and (c) the mutually coupled 
driven oscillators O2 and O3. The symbols PS, IPS, and NPS 
represent the phase synchronization, imperfect phase synchro- 
nization, and unsynchronized states respectively. 



A. Heterogeneous Coupling 

Consider the system of three coupled Rosslcr chaotic 
attractors, 
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We study the different synchronization states in cou- 
pling parameter's space, ei and e 2 . The internal param- 
eters of these oscillators are fixed at ai=0.2, 6i=0.2, and 
ci=5.7 (for driving oscillator Oi ); while a 2 =0. 15, 62=0. 2, 
and c 2 =10 (for mutually coupled oscillators, 2 and O3). 
At these set of parameters all oscillators show chaotic mo- 
tion. The natural frequencies fi turn out to be /1 = 1.079 
and f 2 = f 3 = 1.041 respectively. 

Shown in Fig. [2] are schematic phase diagrams with 
regard to phase synchrony as a function of ei — e 2 . The 
different regimes are characterized through a measure 



based upon the time-dependence of the phase difference 
(see Appendix A) . In order to verify the phase synchrony 
Jj_ 16|we use the phase for Rossler oscillators in Eq. ([1]) 
17l Il8j. namely 4>i = arctan(j/i/xi), i=l, 2, 3. Phase 



synchronization (PS) occurs when the phase difference 
between two interacting oscillators \4>i — 4>j \ saturates [l6| 
(see Fig. [2] for details). In imperfect phase synchrony 
(IPS) regions in Fig. [21 the subsystems are phase locked 
but are subject to occasional slips. The value of J<£i — 4>j\ 
varies during the evolution of the chaotic system [19[ , and 
it changes in a step- wise manner. Each step corresponds 
to a phase synchronized state under a particular phase 
locking condition. The jump between consecutive steps 
occurs in multiples of n In quasiperiodically forced 
systems [2(| the phase differences in the IPS state also 
change in arbitrary multiples of n. Fig. [2] shows phase 
synchronization states between (a) the driving oscillator 
Oi and directly driven oscillator 2 , (b) the driving os- 
cillator Oi and indirectly driven oscillator 3 , and (c) 
mutually coupled response oscillators 2 and 3 . 

Transitions among the phase synchronization states 
when the coupling parameters are varied are depicted in 
Fig. HJa). In one case unsynchronized oscillators (NPS) 
become phase synchronized (PS, region B) via the im- 
perfect phase synchronized state (IPS), and in another 
transition the phase synchronization (PS, region A) tran- 
sitions to the IPS state. As shown in Fig. [2Jb) the indi- 
rectly driven oscillator 3 is initially unsynchronized to 
the driving oscillator Oi , but with the increase of mutual 
interaction between 2 and 3 the oscillator 3 becomes 
phase synchronized to the drive Oi . This phase synchro- 
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FIG. 3: Phase differences between (a) oscillators Oi and O2 
for PS at ei = 1, IPS at ei=2, and NPS at ei=0.15; (b) 
oscillators Oi and O3 for PS at ei=l, IPS at ei=1.2, and 
NPS at ei=0.2; and (b) oscillators O2 and G> 3 for PS at e x =l 
and IPS at £ X = 1.5. The bidirectional coupling parameter is 
fixed at £2 = 1 in all cases. 
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FIG. 4: Phase diagram (schematic) for the coupling parame- 
ters £3 and ei. £3 presents the coupling between O3 and O4. 
The coupling between O2 and O3 is kept fix at £2 = 1. Ci, C2, 
and C3 are separating curves for the pairs O1-O2, O1-O3, and 
O1-O4 respectively. These curves separate PS, IPS, and NPS 
for respective pairs. Details are given in text. 



nization behavior of O3 clearly signifies the transmission 
of drive to the mutually coupled oscillators. Phase syn- 
chronization between the response subsystems O2 and O3 
is seen in Fig. [21(c). By comparing the Figs. HIa),[2Jb), 
and [He) we observe that mutually coupled oscillators O2 
and O3 are phase synchronized when each of the driven 
oscillators phase synchronizes with the drive. 

The above results are further illustrated at different 
points in Fig. [31 Fig. [3J (a) describes oscillators 0\ and 
O2, Fig. |3j(b) describes the behavior of 0\ and O3, while 
the phase synchronization states in between the mutually 
coupled oscillators O2 and O3 is shown in Fig. [3] (c) . The 
bidirectional coupling parameter 62 is kept fixed at £2=1, 
and the different curves are drawn for different values of 
the coupling parameter ei. These curves clearly show the 
behavior of interacting subsystems: the phase difference 
remains bounded for PS whereas it continuously grows 
with time in NPS. 

We further consider a chain of N oscillators with 
nearest-neighbor coupling. Phase synchronization be- 
tween the drive 0\ and responses (Oi, i = 2,3, ... ,N) 
is lost in such systems as we increase the number of 
response oscillators. To determine whether the loss of 
phase synchrony is abrupt or gradual, we first study the 
case of an extended system consisting of a driving os- 
cillator 0\ and three mutually coupled response oscilla- 
tors (Px—^Oi^O^^rO/^). Among three mutually coupled 
subsystems, O2 and O3 are connected by £2 while subsys- 
tems O3 and O4 are connected by the coupling parameter 
£3. Fig. [4] shows the phase synchronization state between 
the subsystem 0\ and driven subsystems O2 , O3 , and O4 . 
To observe the effect of forcing in extended systems, the 



coupling between O2 and O3 is fixed at £2=1, while Fig.0] 
is for varying £3 and ei. 

Our numerical results (see Fig. [4j suggest that phase 
synchronization is not lost abruptly: with the increase of 
coupling between O3 and O4 the single bounded region 
of phase synchronization shrinks between O1-O2, O1-O3 
and the drive and response subsystems turns out to be 
imperfect phase synchronized (IPS) and then phase un- 
synchronized (NPS). The phase synchronization states in 
Fig. [4] are observed for all possible drive-response pairs 
(Oi — O2, 0\ — O3, and 0\ — O4), and in comparison of 
first two pairs we find that the phase space area for phase 
synchronization (PS) is larger in 0\ — O2 then in 0\ — O3 
case. The larger phase space area for PS in 0\ — O2 than 
in O1 — O3 shows that in a chain of mutually coupled oscil- 
lators, the ability of the drive to synchronize the system 
decreases with distance from the drive. 

It should be noted that curves C\ and C2 are bound- 
aries between the phase synchronized (PS) and imper- 
fect phase synchronized (IPS) region for 0\ — O2 and 
0\ — O3 pairs of subsystems respectively (see Appendix 
A). C3 separates the phase unsynchronized state (NPS) 
from the IPS state for the 0\ — O4 pair. The decreasing 
ability of the drive to cause synchrony is further verified 
in Fig. 0] where we see that the oscillators 0\ and O4 are 
largely phase unsynchronized (NPS) in given parameter 
range, though we observe the small appearance of IPS 
state in higher parameter values. The oscillators Oi and 
O4 continue to be in this IPS state for higher values of 
the varying parameters (£3, £1). 
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B. Homogeneous Coupling 

In this scheme two mutually coupled identical chaotic 
systems arc coupled unidircctionally to another (but non- 
identical) chaotic oscillator, unidircctionally through the 
same variable (here y). See Fig. [TJb). The equations 
of motion for three coupled Rosslcr oscillators in this 
scheme are 

x\ = -yi - zi 

ill = xi + aiyi 

zi = h + z±(xi - ci) 

%2 = -V2-Z2 

in = x 2 + a 2 y2 + £1(2/1 - 2/2) + £2(2/3 - 2/2) 

%2 = h + z 2 (x 2 - c 2 ) 

X3 = -2/3 - 23 

2/3 = £3 + a 2 y 3 + £2(2/2 - 2/3) 

Z3 = b 2 + 23(2:3 - c 2 ). (2) 

A schematic phase diagrams as a function of the pa- 
rameters ei and e 2 is shown in Fig. [SJ Fig. [5] (a) shows 
the phase synchronization states between the driving os- 
cillator 0\ and directly driven oscillator O2, while Fig. [5] 
(b) shows the phase synchronization states between the 
driving oscillator 0\ and indirectly driven oscillator O3. 
The phase synchronization states between the mutually 
coupled response subsystems O2 and O3 are shown in 

Fig. m (c). 

The occurrence of phase synchronization states be- 
tween 0\ and O2 or O3 is a consequence of the trans- 
mission of forcing through mutually coupled chaotic os- 
cillators. In order to compare this effect with that of 
the heterogeneous case (previous subsection) we consider 
the larger number of oscillators in a chain. As shown 
in Fig. |5] by increasing the number of oscillators to 4 
(Oi — > 2 f-» O3 f-> O4) with the coupling between 2 
and O3 fixed at e 2 — 1 we find that the external drive 
transmission decreases with distance but the phase syn- 
chronization between the drive 0\ and responses 2 , O3, 
and O4 continues to occur so long as the coupling is het- 
erogeneous. These results hold for even larger numbers 
of mutually coupled oscillators. A second difference be- 
tween hetro- and homogeneous schemes is that there are 
multiple regimes of PS in the former case (see Fig. [2]) 
but only a single regime for latter. 

III. FORCING THROUGH MEDIATING 
SYSTEMS 

Natural systems are often modulated indirectly. Con- 
sider a drive-response pair mediated by a number of mu- 
tually coupled subsystems. In order to clearly under- 
stand the present case of external forcing via intermediate 
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FIG. 5: Schematic phase diagram for the synchronization 
states as a function of the coupling parameters ei and £2 in 
the homogeneous coupling scheme between (a) oscillator 0\ 
and the directly driven oscillator O2 , (b) 0\ and the indirectly 
driven oscillator O3, and (c) mutually coupled oscillators O2 
and O3. The symbols PS, IPS, and NPS have the same mean- 
ing as in Fig. [2] 



subsystems, the model system is crafted from both het- 
erogeneous and homogeneous coupling schemes. (Recall 
that the difference between these two coupling schemes 
is as follows: homogeneous coupling has both the drive 
and response subsystems coupled by the same variable, 
while heterogeneous coupling uses different variables.) 

The heterogeneous coupling is shown in Fig. [7] (a) 
where oscillator 0\ drives 2 and O3 which then drive os- 
cillator O4. If £3 = this reduces to the case of Fig. [TJa) 
but for finite £3 different phase synchronization states 
(PS, IPS, and NPS) are observed between oscillators 0\ 
and O4. Results are shown in Fig. [7| (b) for a repre- 
sentative value of £2 = 1. As the number of mutually 
coupled oscillators between the concerned drive-response 
pairs (e.g. 0\ — O4) is increased, the PS regime is lost. 
This result though not presented here, is in consonance 
with the result of Fig. [1] for the heterogeneous coupling 
case. 

With homogeneous coupling, as shown in Fig. [7] (c), 
phase synchrony between the terminal oscillators 0\ and 
O4 is also achieved. Results are shown in Fig. [JJ(d) and it 
can be seen that for a given value of £3 oscillators 0\ and 
O4 become phase synchronized beyond a critical value 
of coupling strength ei . The non-monotonic form of the 
border between IPS and PS in Fig. [7J (d) is possibly due 
to two indirect forcing: 0\ is forcing 2 and O3, while 
the combined oscillators Oi,02, and O3 are forcing O4. 
This PS regime persists even in case of increased number 
of mutually coupled intermediate oscillators. 

Heterogeneous and homogeneous coupling schemes 
thus clearly exhibit distinct dynamics. In the former, 
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FIG. 6: Schematic phase diagram for the extended case, show- 
ing synchronization states as a function of coupling parame- 
ters £3 and ei. £3 is the coupling between O3 and O4, £1 
couples the drive to the responses. Phase synchronization 
states between (a) drive 0\ and the directly driven oscillator 
O2 (b) 0\ and the indirectly driven oscillator O3, and (c) 0\ 
and the indirectly driven oscillator O4. In all cases, £2 = 1. 



synchronization regimes are truncated while in the later 
case synchrony persists even when the number of mu- 
tually coupled oscillators is increased. The conjugate 
variables employed in heterogeneous couplin g ap pear to 
provide an effective time-delayed interaction [12, [2IJ . 

By computing the conditional Lyapunov exponents 
[2^ | as a function of the external coupling strength, we 
find that in general, phase synchronization is not ob- 
tained when the drive couples to z, so in the present 
studies, we drive the variable y (see Eqns. ([T]) and ((2J). 
Results are similar when drive is coupled with the x vari- 
able as well. 



IV. SUMMARY 

In the present work we have studied the effect of an 
external drive on symmetrically and diffusively coupled 
response oscillators. Two different coupling schemes have 
been examined; these coupling schemes were based on the 
use of the two coupling parameters either on the same 
variables (homogeneous coupling) of interacting oscilla- 
tors or on different variables (heterogeneous coupling). 
Numerous combinations are possible, and this study is 
an attempt to consider some representative examples. 

The relaying of the drive to both the directly and the 
indirectly forced subsystems is seen through the occur- 
rence of phase synchronization. For heterogeneous cou- 
pling phase synchronization between drive 0\ and re- 
sponses Oi, i = 2, 3, ...,N is lost in extended systems, 
but in a gradual manner. Thus the transmission effect of 
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FIG. 7: Heterogeneous coupling: (a) the system of Fig. [T] (a) 
drives a dissimilar oscillator O4 through the response oscilla- 
tor O3, (b) a schematic phase diagram as a function of £3 and 
£1, between oscillators Oi and O4 for fixed £2 = 1. Homoge- 
neous coupling: (c) all oscillators are connected through the 
same variable, here y, and (d) phase synchronization states 
between oscillators 0\ and O4 for fixed £2 = 1. 



external drive decreases sequentially in an array of mu- 
tually coupled response oscillators. In case of homoge- 
neous coupling the external drive is transmitted to each 
of the mutually coupled response subsystems. The effect 
of indirect external forcing has been also observed for a 
typical model with mediating oscillators. 

The occurrence of phase synchrony is of particular sig- 
nificance due to potential applications in diverse fields 



G 



such as biological systems flij . networks of coupled dy- 
namical systems [Iff among others. The present results 
indicate that the transmission of drive to an array of 
mutually coupled response oscillators depends strongly 
on the coupling pattern, and it is very likely that the 
topology of coupling will also play a significant role. A 
study of different coupling motifs in order to determine 
the manner in which the drive is transmitted in more 
complex networks is therefore presently under way [23j . 
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Appendix A: Phase Synchronization 

The manner in which the phase difference varies with 
time distinguishes the states of PS, IPS, and NPS. When 
two interacting systems are phase synchronized, the 
phase difference between the systems remain bounded Q 
while the phase difference between the interacting sub- 
systems continuously grows when the systems are phase 
unsynchronized (NPS). In case of imperfect phase syn- 
chronization the phase difference grows in multiples of 
7r [l9|, [2(| ■ This variation in phase difference is used to 



deduce a measure for the identification of different phase 
synchronization states. The averaged phase difference 
increment $ is the time average of the derivative of the 
instantaneous phase difference A0 t = | </> 2 — <f>\ | at time 
i, namely 

$ = (A& +1 -A&). (Al) 

1. If the variation of phase difference remains bounded 
between and n, $ should be close to zero in case 
of phase synchronization. 

2. For IPS the phase difference between the two in- 
teracting subsystems grows in multiples of 7r; this 
results in large values of <F 

3. When the interacting subsystems are out of phase 
synchrony, the phase difference increases with time 
but <f> in considerably lower than it is in the case of 
IPS. 

Fig. [8] shows the averaged phase difference increment 
$ as a function of the parameter e\ for the system in 
FigfT] (a) with fixed e 2 = 1. The inset shows the regions 
of PS, IPS, and NPS. This measure is qualitative in the 
sense that sharp boundaries between these states (Figs. 
[H HI El ElEJb) and (d)) cannot be drawn in a quantitative 
manner. We use thresholds for $ as 10~ 5 , 0.0001, and 0.1 
to assign the boundary for the transition PS to IPS/NPS, 
NPS to IPS, and NPS to IPS respectively. 



[1] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 64, 821 
(1990). 

[2] A. S. Pikovsky, M. G. Rosenblum, and J. Kurths, Syn- 
chronization: A Universal Concept in Nonlinear Sciences 
(Cambridge University Press, Cambridge, U.K., 2001). 

[3] K. Kaneko, Theory and applications of coupled map lat- 
tices (John Wiley and Sons, NY, 1993); E. Ott, Chaos in 
Dynamical Systems (Cambridge University Press, Cam- 
bridge, 1993); A. Prasad, L. D. Iasemmidis, S. Sabesan, 
and K. Tsakalis, Pramana-J. Phys., 64, 513 (2005). 

[4] H. Fujisaka and T. Yamada, Prog. Theor. Phys. 64, 821 
(1983). 

[5] A. S. Pikovsky, Z. Phys. B 55, 149 (1984). 

[6] N. F. Rulkov, M. M. Suschik, L. S. Tsimring, and H. 

D. I. Abarbanel, Phys. Rev. E 51, 980 (1995); H. D. I. 

Abarbanel, N. F. Rulkov, and M. M. Suschik, Phys. Rev. 

E 53, 4528 (1996). 
[7] M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, Phys. 

Rev. Lett. 76, 1804 (1996); M. G. Rosenblum, A. S. 

Pikovsky, and J. Kurths, Phys. Rev. Lett. 78, 4193 

(1997). 

[8] L. Kocarev, and U. Parliz, Phys. Rev. Lett. 76, 1816 
(1996). 

[9] A. Prasad, Chaos, Solitons & Fractals 43, 42 (2010). 
[10] Y. Kuramoto Chemical Oscillations, Waves, and Tur- 
bulence (Springer, New York, 1984); R. E. Mirollo, and 



S. H. Strogatz, SIAM J. Appl. Math. 50, 1645 (1990); 
N. Kopell, G. B. Ermentrout, M. A. Whittington, and 
R. D. Traub, Proc. Natl. Acad. Sci. U. S. A. 97, 1867 
(2000); T. J. Lewis and J. Rinzel, J. Comput. Neurosci. 
14, 283 (2003); E. M. Izhikevich, Dynamical Systems in 
Neuroscience: The Geometry of Excitability and Burst- 
ing (MIT Press, Cambridge, MA, 2007); A. Sherman, J. 
Rinzel, and J. Keizer, Biophys. J. 54, 411 (1988); M. G. 
Pedersen, R. Bertram, and A. Sherman, Biophys. J. 89, 
107 (2005). 

[11] A. Prasad, M. Dhamala, B. M. Adhikari. and R. Ra- 
maswamy, Phys. Rev. E 81, 027201 (2010). 

[12] R. Karnatak, R. Ramaswamy, and A. Prasad, Phys. Rev. 
E 76, 035201 (2007). 

[13] R. Karnatak, R. Ramaswamy, and A. Prasad, Chaos, 19, 
033143 (2009). 

[14] J. D. Murray, Mathematical Biology (Springer- Verlag, 

New York, 1993). 
[15] M. Barahona and L. M. Pecora, Phys. Rev. Lett. 89, 

054101 (2002); H. Hong, M. Y. Choi, and B. J. Kim, 

Phys. Rev. E 65, 026139 (2002); C. Li, W. Sun, and J. 

Kurths, Phys. Rev. E 76, 046204 (2007). 
[16] M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, IEEE 

Trans. Circ. Sys. 44, 874 (1997). 
[17] A. S. Pikovsky, M. G. Rosenblum, and J. Kurths, Euro- 

phys. Letts. 34, 165 (1996). 



7 



0.2 



© 0.1 




[18] A. Goryachev and R. Kapral, Phys. Rev. Lett. 76, 1619 
(1996). 

[19] E.-H. Park, M. A. Zaks, and J. Kurths, Phys. Rev. E 

60, 6627 (1999). 
[20] M. Agrawal, A. Prasad, and R. Ramaswamy, Phys. Rev. 

E 81, 026202 (2010). 
[21] N. H. Packard, J. P. Crutchfield, J. D. Farmer, and R. S. 

Shaw, Phys. Rev. Lett. 45, 712 (1980). 
[22] U. Parlitz, L. Junge, W. Lauterborn, and L. Kocarev, 

Phys. Rev. E 54, 2115 (1996). 
[23] M. Agrawal, in preparation. 



FIG. 8: Phase synchronization measure for heterogeneous 
coupling with £2 = 1, varying ei. The inset shows the or- 
der of the averaged phase difference increment in the three 
different states. 



